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Abstract 

The calculation of the radiative energy loss encountered by a fast charged particle which under- 
goes multiple scattering is being investigated. A detailed derivation of the Landau- Pomeranchuk- 
Migdal coherent effect in QED is given, focusing on the specific feature of the Coulomb interac- 
tion. As a result the radiation intensity per unit length in the coherent regime is shown to be 
proportional to ^/LJ for a photon energy u times a logarithmic enhancement which is determined 
exactly. 
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1 Introduction 



The radiative energy loss encountered by a charged fast particle which undergoes successive 
scatterings in a medium has recently been the object of renewed interest. This is indeed of 
special importance in order to understand the energy loss mechanism for quarks and gluons, 
propagating in nuclear matter and in particular through the quark-gluon plasma. 

Our first step ^ has been to reformulate the QED treatment given by Landau, Pomeranchuk 
and Migdal (LPM) [|,|,|]. The experimental observation of the LPM effect at SLAG [|] has 
recently triggered further studies with special emphasis on finite length effects p. 

In the present paper we focus on a detailed derivation of the radiative spectrum per unit 
length, assuming a very large number of scatterers. The model used which has been recently 
elaborated by Gyulassy and Wang ^ depicts the multiple scattering of a fast electron in the 
medium as due to static scattering centres with Debye screened Coulomb potentials. The as- 
sumption that the mean free path A of the projectile is much larger than the screening radius 
A ^ allows one to treat successive scatterings as independent. Coupled to the soft photon 
approximation this then leads to an eikonal picture of classical propagation. 

As stated in [|I|, the specific case of the Coulomb potential which is not screened at short 
distances requires a special treatment which was not available in the literature. Such a treatment 
is essential for QCD. We present here for QED a complete detailed derivation, correcting an 
unjustified step made in [^. 

For the specific case of the Coulomb interaction, the radiation intensity per unit length in 
the coherent regime is proportional to a/cJ times a logarithmic enhancement. For potentials less 
singular than the Coulomb one at short distance, where the random walk picture is applicable, 
we confirm the original LPM result. 

The outline of the paper is as follows : 
In section |^, we show in details how the model of static Coulomb centres describes multiple 
scattering. The radiation spectrum induced by multiple scattering is worked out in the soft 
photon approximation in section ^ and studied in section ^. In this last section, we first obtain 
the shape of the spectrum on heuristic grounds and then carefully derive the spectrum in the 
coherent regime by solving an exact differential equation to leading logarithmic accuracy. Section 
^ is the conclusion. 

2 Model for multiple scattering 

In order to work out the radiation intensity properly normalized we need first to derive the 
multiple scattering cross section of a fast electron propagating in a medium. We do this using 
the model of Gyulassy and Wang 0. The main feature of this model consists in assuming that 
scattering centres are static. Thus the coUisional energy loss of the charged particle vanishes, 
and the total energy loss will be purely radiative. The scattering centre located at Xj creates 
screened Coulomb potential 





with Fourier transform 



e 



-iq-Xi 
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Here /i is the Debye mass in the medium and we denote the momentum transferred to the 
electron during scattering on the centre Xj. 

To be able to treat successive elastic scatterings as independent, centres have to be well 
separated. This means that the average distance between two successive interactions (the electron 
mean free path A), is large compared to the range of the potential, 

A>/i-^ (2.3) 

This leads to important simplifications in the study of the radiative energy loss. 

We also concentrate on the high energy limit, where the transferred momenta qi become 
transverse. Moreover, their characteristic values are of order /i, and we have 

e. = ^ ~ = I « 1 , (2.4) 

where 0s is the typical electron scattering angle. 

Let us consider the probability amplitude to scatter on N static centres, and call the incoming 
and outgoing electron 4-momenta po cind ptv, respectively. In the high energy limit where electron 
mass and spin effects are irrelevant, the S"- matrix element corresponding to multiple elastic 
scattering reads g 

Sscatt QC 6{p% - Pq)Y, n 
o- i=l 

where after the i*^ momentum transfer g^, the electron momentum is 

Pi=Pi-i + qi ; i = l,...N, (2.6a) 

and 



d^Pi 
pf + ir] 



n 

i=l 



(2.5) 



pO = E ; i = 0,...N -1 . (2.6b) 

The resulting cross section may be worked out in the eikonal limit of small angle scattering (see 
e.g. 0). Let us indicate the main steps of the derivation in the well known Glauber picture. 
Sscatt contains the phase factor exp{iipscatt) with 

N N N~l 

i=l 1=1 i=l (2-7) 

+ {P0\\ -PN\\ Za{N)) , 

where Zo-(i) is the longitudinal component of Xo-(i) and we sum over all permutations a of {1, . . . A^}. 

Sscatt may be calculated by integrating successively over the longitudinal momenta Pi\\, for 
i = 1, . . . N — 1. Let us start with / dpi\\. 

When {Zo-(2) — Za{i)) > 0, the integration over pi\\ is performed by closing the integration 
contour in the upper half of the complex pi\\ plane (see ( p.7| )), and taking the pole = E — 
pIJ2E + irj. 
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When {z„(^2) ~ ^cr(i)) < 0, we have to close the contour in the lower half-plane. We take the 
pole pi\\ = —E + pIj_/2E — ir], but this yields a residue which is suppressed by a factor of order 
{p1j_ + ij?)/E'^ when compared to the first case, where we have chosen pq± = 0. 

Finally, whatever the sign of (;Zo-(2) — is, the poles in pi\\ corresponding to the "Coulomb 
denominator" (gl^ + /i^) lead to residues involving a suppression factor ~ exp{— /i|zo-{2) — -2o-(i)|}, 
and may be neglected when A ^ This may be repeated for all Pi\\ integrations. 

We deduce from the above : 

• The high energy electron scatters successively on ordered scattering centres (in the z di- 
rection) : Zo-[N) > ^fr(iv-i) > ••• > ^o-(i)- In (p.5|), only the identity permutation has to 
be kept and backward scattering may be neg lectedQ. 



As a consequence of ( p.3|) , the internal electron momenta are on-shell, which corresponds 
to independent elastic scatterings on static centres: 

Pill^i?-^, 2 = 1,2,. ..iV-1; Po\\=E. (2.8) 



The phase 



takes the form 



Vscatt 



2E 



N N-1 p2 

i=l i=l 



(2.9) 



Using Qij^ as integration variables we obtain, up to an irrelevant phase and the energy conservation 
(5-factor, 



Sscatt i W_ 



N 



i=l 



-iqi±-Xi 



7V-1 

n 

i=l 



N 



(2.10) 



The scattering amplitude squared reads 



\Mscatt\'^ OC 



N 

n 

i=l 



N-l 

n 



^1=1 



2E (zi+i-Zi) 



(2.11) 



The next step is to average over the coordinates Xj. Since the centres are assumed to be uniformly 
distributed, averaging over Xi± leads to q[j_ = qi±, p[^ = Pi±, which results in 



OC / \[ 



d^p 



•Nl. 



N 



=1 {Qil + ^^ 



2\2 



\ PNi 



(2.12) 



\i=l 



Here we have dropped an overall factor including the transverse surface 5^(0). Integrating over 
Pni. we finally arrive at 



A similar expression may be found in P,pi|. 

^In other words the permutations different from identity imply the occurence of vacuum creation of e^e~ 
pairs, which process is power-suppressed when E oo. 



4 



3 Radiation spectrum induced by multiple scattering in 
the soft photon approximation 

Multiple scattering induces radiation. The calculation of the bremsstrahlung probability is done 
in the soft photon approximation 

oj <^ E . (3.1) 

This allows to factorize the radiation amplitude as the product of the multiple scattering am- 
plitude times the photon emission amplitude. The ratio between the radiation cross section and 
the scattering cross section yields the radiation probability. 

Let us consider the probability amplitude for emitting a photon of 4-momentum k = {uj, k) 
off the electron between centres Xj and Xj+i. The corresponding 5*- matrix element reads 



i=l 



pf + ie 



.3 - ^ ■ Pj 
k-pj 

N~l 

n 

i=j+i 



{pj — ky + ir] p'j + it] 
d^Pi 



{pi - ky + iri 



N 

n 

i=l 



(3.2) 



with the momentum variables pi defined in ( |2.6|) . The two physical photon polarization states e 
are chosen in the form 



e-k = 



e± ■ k± e±- kj 

u + k\\ 2uj 



(3.3) 



In the same way as in the previous section, the phase (frad is given by the r.h.s. of ( |2.7| ) and the 
integrations over pm are done by closing the contour in the upper half-plane, 



Pi\\ 



2E 



{P^ - kf 



Piu — k\\ E — UJ 



{P^ - k)l 

2{E-uj) 



(3.4a) 
(3.4b) 



As far as the intermediate state j is concerned, we may have either {pj — k^ = or p| = 0. Since 
the accompanying radiation factor is invariant, 

e-pj _ e- {pj - k) _ 



k ■ Pj k ■ {pj — k) 



e-k = 0, k^ 



(3.5) 



it may be always expressed in terms of the real electron momentum. For example, taking pj = 
we obtain 



6-pj 

k-pj 



£o{E + Pj\\) -e±-pj_ 

2 



Ee^- 



k± Pj± 
u E 



= E e± - Uj , 



Euj I k I 



Pj± 
E 



Ecu n 



(3.6a) 
(3.6b) 
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Here we introduced the dimensionless transverse vector Uj which, in the quasi-colhnear kinematics 
under interest, measures the angle between the photon and the electron after scattering number 
h 



u E 

For another pole, {pj — = 0, one has to substitute in ( pl6| 



(3.7) 



^ _ Pj- 



E 



111 . 
— ' 



(3.8) 



which difference, however, may be neglected in the soft approximation ( p.l| 
Thus, we can write the ^-matrix element (13. 21) as 



N-l N 



i=l 1=1 

where the elementary soft radiation factor is 



k-pj 



e-pj 



k-pj 



e-pj 



u 



Pi=0, 



(3.9) 



'J 

— — £± ■ ^ , 



(3.10) 



with Uj defined in (|3^). 

Now we must express the phases (prad- 



^rad 



N 
i=l 



N 



N 



i=l 

N-l 

PO\\Zl -Pn\\Zn + Pi\\i^i+1 - ^i) 



(3.11) 



i=l 



According to (|3.4|), the value of the longitudinal component Pi\\ depends on whether = or 
{pi — kY = 0. Expanding ( |3.4b| ) m uj / E and neglecting 0{uj'^ / E) for the second case, we derive 



Pi 



E 



Pii^ 
2E 



UJ 9 



for {p, -kY = 0, 



(3.12) 



Using (|3.4a|) and (|3.12|) and invoking ( p.9| ) for the scattering phase, we represent the phase of 



(3.13) 



the first term in (|3.9| ), corresponding to {pj — k) = 0, as 

'Prad = 'Pscatt + , 
N-l 



UJ 



2 - 



1=] 



Similarly, for = 0, we use ( p.4aD for i = 1, . . . j and ( p. 12]) for z = j + 1, . . . — 1 to obtain the 
expression for iprad, which differs from ( p.l3|) simply by substituting j + 1 for j, that is $j —>■ $j+i 
in (|3l3|) . 
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We stress that the phase difference depends only on the longitudinal separations. 

Using the picture of the electron moving along a classical trajectory between scattering cen- 
tres, Xj+i — Xi = Vi (tj+i — ti) with Vi the electron velocity, one may write the phase in Lorentz- 
invariant form as 

~ ■ {xj - xn) ■ (3.14) 

After summing S^^^ over j and adding the radiation off initial and final lines we get the total 
radiation S'-matrix element 



N 



(3.15a) 



where 



7V-1 



k-po 



N 

E 



''P' ^■^^•-Mei*.. (3.15b) 



Thus, the integrand of Srad is proportional to that of Sscatt in ( p.lO|) . As seen from (|3.13|) , 
the proportionality factor eTZ does not depend on the transverse coordinates Xi±. Therefore, 
after squaring ( 3.15 ) to get iM^adp, we may follow the same averaging procedure over Xi± as for 

I Mscatt I 



The result is 



N 



(3.16) 



pol 



Integrating over pn± gives 



darnd oc 



.=1 {%l + 1^ 



pol 



(27r)32a; 



(3.17) 



Finally, using ( p.lO| ) and summing over the photon polarizations leads to the radiation intensity 



UJ- 



dl a 

duj TT^ 



dVL 



N 



1=1 



(3.18a) 



with dVt representing the integration over the photon direction. Here we have defined the emission 
current 



A, 



(3.18b) 



The brackets denote the averaging over transverse momenta and longitudinal coordinates using 
normalized probability distributions. 



)) 



Af-l 

n 



d^o 



exp 



A, 



^ ■ n 



N 



=1 T^iQil + 



2\2 



(3.19) 
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where the exponential factor reflects the survival probability of the electron over the distance 
= Z£+i — ze- 

In a general case of scattering potential other than Coulomb, the averaging over transverse 
momenta has the form 



N 

I[d%,V{ql) {...)■ £qV{q') = l 

i=l 



(3.20) 



with V{q^) the normalized cross section for elastic electron scattering in the medium. 

For the sake of completeness and also in view of the generalization to nonabelian radiation 
processes, it is useful to mention the old-fashioned perturbation theory proof of ( p.l8| ). 

Let us consider the radiation amplitude induced by transfers of fixed momenta at centres 
Xi and times ti. In time-ordered perturbation theory, in the hmit E ^ oo, is obtained by 
integrating the phase factor over the emission time Tj, lying between tj and tj+i, 



-le 



2g • P3 
2E 



-j+i 



scatt 



(3.21) 



We denote by 6ip the phase difference between M^^^ and Mg^att- For the interaction times t. 
i > j, the associated phase is ti{\pi \ — \pi-i\) — for Mgcatt, but becomes ti{\pi — k\ — Ipi-i — k 
ti k ■ {vi-i — Vi) for M^^^. Thus 
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6(p 



N 

i=j+i 



[Vi-i - Vi 



Using 
and 

we easily get 



3^1+1 ~ Xi — Vi (t 



\Pi - k\ - \pi\ + U C:^ u - k ■ Vi 



k ■ Pi 



E 



MLd 



k ■ pj ^ 



scatt 



(3.22) 

(3.23) 
(3.24) 

(3.25) 



Integrating over tj in ( |3.21| ) gives two contributions, proportional to e*'^'^^ and e^^'^^+'^, which 
may be respectively included in radiation amplitudes induced by the j*'* and (j + 1)*'^ scatterings. 
Finally, the radiation amplitude induced by scatterings is found to be identical to the previous 
result (CT) . 



4 Explicit calculation of the radiation intensity 



Here we present the analytic calculation of the energy spectrum (|3.18|) in the soft photon ap- 
proximation. The spectrum (|3.18| ) may be written in two equivalent forms as follows: 

N N N 



dl 

djjj 



vr 



I d^ (2 ReE E e**'" + E 

\ i=\ j=i+l i=l 

/ N N 

dn ( 2 ReE E ^i-^j (e'*'" - 1) + 



i=l j=i+l 



N 



EA 

i=l 



(4.1a) 
(4.1b) 
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with the relative phase 



i-i 

$ .. = - = ^ ^ uj{ze+i - ze) . (4.2) 



4.1 Heuristic derivation of the LPM spectrum 

The radiation pattern depends crucially on the phases. If $jj are large, the rapid oscillations 
in the first term of ([4.1a|) wash this contribution away, and we have the Bethe-Heitler regime 
of independent radiation. Dividing by the size of the medium, A^A, for the differential energy 
spectrum per unit length (the radiation density) we obtain 



dl 1 a /• _ 1 



dudz A TT^ 7 



N 



\i=\ 



a f d^ui 



a 



d^u I \ 2a I 62 

/ \ ^, / In 



(4.3) 



TtA J TT \y2(^;^_Q^^2 I 7rA\ U, 



2 

rain i 



with -C 0^ the value of the emission angle below which the phase between neighbouring 
centres becomes small: 

= \ujul{zi+i - z,) ~ I^Am- < 1 ^ m- < ; = i^A. (4.4) 

Small emission angles u"^ < or (u — Qs)^ < u^in correspond to photon formation time larger 
than the distance to the neighbouring centres. 



2uj ^ . , 

tform = ^ = ^ > (4.5) 



SO that destructive interference with the radiation due to the nearest neighbour screens the 
logarithmic divergence. Introducing the characteristic dimensionless parameter k, representing 
the typical phase difference between neighbouring centres. 



1^ = ?r 

2 

we may obtain (4.3) as 



(4.6) 



The Bethe-Heitler limit corresponds to k > 1. 

In the opposite limit, k -C 1, the phases between neighbouring centres are vanishingly small, 
so that a group of u centres radiates coherently. To estimate the coherence number v we 
look for the separation between two centres, such that an accumulated phase becomes of order 
unity. In a first approximation one may assume that the trajectory of the electron is a random 
walk in the transverse momentum space, 

^ u\ + (en . (4.8) 
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The phase ^i^i+u is estimated then as 



uul + y(u - 1) (Q 



el 



el 



(4.9) 



As will become apparent later from ( [4.27a| ), the second term here is v times bigger than the first 
one, which gives 



K (92) 



\2 ei 



(4.10) 



If the coherence length exceeds the size of the medium, z/A > L, all the phases are small, and 
the first term in ( [4.1b| ) vanishes. In the last term the emission currents ( |3.18b| ) add up into the 
expression 

Qtot 



N ^ ^ 

^ X UN Uq 

Z^Ai = — 2 ' UN = Uo 

i=l Uq 



(4.11) 



which corresponds to the radiation induced by a single scattering with the momentum transfer 



N 



(l±tot = '^Qi- It is usually called the factorization limit. For a given qtot the radiation 



IS 



1=1 



independent of the size and the properties of the medium and reads 



duo 



TX 



(4.12) 



Here we have introduced a finite electron mass m to regularize the collinear divergence. It is 
the only place where the electron mass enters in the high-energy limit. Otherwise, the induced 
radiation is coUinear-safe. When the size of the medium is large enough to embody a few 
coherence lengths, L ^ z/A, the LPM-suppression of the Bethe-Heitler spectrum takes place. In 
order to quantify it, it suffices to "slice" the medium and substitute the number of "effective 
radiators" N/ v for A^. For the radiation density we write 



a; 



dl \ I dl ^ 

~ — ijj 

dwdz V \ dujdz 



(BH) 



(4.13) 



In the standard Bethe-Heitler spectrum off a point-like source (|4.3| ) the logarithmic enhance- 
ment factor appears due to emission angles much smaller than the scattering angle. The effective 
centre, however, radiates at typical angles such that the formation time 2{uju'^)~^ is of the 
order of the length of the radiator z/A. This implies an emission angle of the order of the 
accumulated diffusion angle, 

^ (4.14) 



Under these conditions the logarithmic enhancement is absent, and 13|) becomes 



dl 



1 2a 



UJ 



duodz V ttA 



C = C 



a 



ttA \ 



2k 



c 



a 



TX 



AE2 



C 



a 



CO — ^ — 



TT 



(4.15) 



10 



with C a constant of order 1. This reproduces the well-known Landau- Pomeranchuk-Migdal 
result (C = 1). Strictly speaking, the above derivation based on the random walk picture applies 
only to the cases when the mean squared momentum transfer in a single scattering is well defined. 
However, this is not true for Coulomb scattering where the integral determining {q\) formally 
diverges logarithmically: 

d)-/7#^-.= -.'/7#^-.^ln%. (4,6) 



Replacing the upper limit of the angular integral by the characteristic angle of the problem v? 
from ( 4.14| ), we obtain 



> - I In-, (4.17) 



2 K 



where we have neglected the log-factor under logarithm. Substituting this ratio into ( [4.15|) , for 
the Coulomb case we derive 

( dl \ a I ^ 1 ^ 

This heuristic estimate coincides with the true answer derived below. 



The origin of the extra logarithmic enhancement in (|4.18|) is readily understood. The singu- 
larity of the Coulomb potential at small distances corresponds to a long tail in the transverse 
momentum distribution. This enriches the contribution from "large jumps" , when the momen- 
tum transfer exceeds the inverse Debye radius, q]_ ^ /i^. Effective "random walk" steps become 
(logarithmically) larger, which reduces by the Vln factor the coherence number v and, thus, 
the LPM suppression. 

Hereafter we shall concentrate on the coherent LPM regime 

I <^ N =^ -^<K<1. (4.19a) 

In terms of the photon energy the coherent region is limited by 

< — < ^ = 7^ • (4.19b) 



cr 



4.2 General expression for the induced radiation spectrum 

We turn to the original expression for the radiation spectrum ( [4.11 ). As we shall see below, a 
finite (though large) number of scatterings (j — i ~ v) is essential, so that for N ^ u a. given term 
of the double sum clearly depends only on the relative position of the two centres, j — i = n + 1, 
n > 0. Moreover, the internal sum over j in ( [4.1b| ) converges and, therefore, does not depend on 
i in the same approximation. The sum over i then gives the total number of scatterings N. 
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Dividing by the size of the medium, NX, and neglecting the factorization contribution in the 
N —>■ oo hmit, we obtain the following expression for the differential energy spectrum per unit 
length: 



UJ 



dl 



dujdz TiX J 71 



2 Re ^ Ji ■ Jn+2 

\ n=0 



n+1 



exp iiKy^ U, 



2 ^e+i ~ 



l + O(-))(4.20) 



Here we have introduced new variables Ui to represent ui in units of the typical scattering angle 
9, = fi/E : 

= . ; dQ = el d^tj , (4.21) 
and rescaled the emission currents correspondingly: 



It is also convenient to express the transferred momenta in units of /x, 



Qt 



1^ 



Ue — Ui-i — Qi 



(4.22) 



(4.23) 



so that 



Qe — Ug — Ue-i 



1 . 



Performing in (|4.2CI|) the averaging over longitudinal separations with use of (|3.19| ), we obtain 



UJ- 



dl 

dujdz 



71 A J n , nJ n 1 



oo „ n+2 

E/n 



n+1 



n ^(^m) - 1 



.m=l 



where 



2\-l 



(4.24a) 



(4.24b) 



This holds in general for arbitrary interactions. For the particular case of Coulomb scattering, 
V is given by 



— * — » 

In ( [4.24| ), the dependence on Q\ and (5n+2 is contained only in the product of currents 



(4.25) 



Un+l — Q 



n+2 



n+2 



Un+1 

UI {Ul+QlYl UUn+l-Qn+2? Ul+1 



(4.26) 



Keeping f/i, U2, ■ ■ - Un+i fixed, we integrate first over Qi (which is equivalent to integrating 
over Uo, the direction of the incoming electron with respect to the photon) and Qn+2 (outgoing 
electron). To this end we define 



fo{U,) ^ I d'Q,V{Qi) Ji = ^-^J dQi V{Qi) &{Qi - Uf) = tt^ J^^ dQ' V{Q') . (4.27a) 
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In the same way we have 



d^Qn+2V{.Qli+2)Jn+2 



TC 



Un+1 

Un+1 JUI^^ 



(4.27b) 



The spectrum can be rewritten as 

dl 2a 



LU- 



— Re 



duj dz rrX 



d^Ui 



TX 



foiUi) ■ fVi 



K=0 



(4.28) 



where f{Ui) is given by 

oo n+1 p ^ 

d'QeViQ',)ijiUf) 
The K-dependence of / comes from the ^-dependence of ip (see ( [4.24b| )). 



oo n+1 
n=l e=2 



n+1 J 



(4.29) 



Master equation. By expanding the sum over n, it is straightforward to see that f{U) satisfies 
the following integral equation 

(4.30) 



(1 - tKU')f{U) = fo{U) + J d'Q V{Q') fill - Q) 
For K = it has a trivial potential-independent solution 

Indeed, substituting (|4.31|) under the integral, we have 



U 
IP 



(4.31) 



d'Q Vm -^S- = vr^ dQ' V{Q') 
^ (f/ _ g )2 Jo 

Taken together with ( ^.27| ) for /o, this gives for the r.h.s. 

U 



IT 



[/2 



{/•CO rU^ 
/ dQ^ ViQ^) + / dQ^ F(g2 
Jm Jo 



u_ 
IP 



■ 1 



(due to the normalization of the scattering cross section ( p.20|) ), which is identical to the l.h.s. 
(«: = 0). 



In order to solve the equation ( |4.3CI| ) for k 7^ we adopt a method different from the derivation 
advocated in fll], which relied on an incorrect approximation. Going to the impact parameter 
space by defining the Fourier transform 



f{B) 



d^U e-*^-^/(f7) , 
d^Q e"'^-*^ V{Q^) 



(4.32a) 
(4.32b) 
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we first derive the -B-space image of tlie function ( 4.27 ): 

B 



MB) = -2m^il-ViB')). (4.33) 

In these terms ( [4.30| ) converts into the differential equation 

(l + mV^)/(S) = foiB) + ViB') fiB). (4.34) 
Now we introduce the scalar function g 

7iB) = ^ ~giB') (4.35a) 



to obtain 



V|?(5) = 4 B g"{B') , g" = . (4.35b) 

Finally, defining 

h{B^) = g{B^) + 2TTi = B ■ f{B) + 2m , (4.36) 
we represent (|4.34|) in the form 



AiK h"{B^) + - — h{B') = . (4.37a) 
B 

This is a linear second order differential equation. The corresponding boundary conditions are 



/i(0) = 27rz, (4.37b) 

h{oo) = 0. (4.37c) 

Indeed, ( [4.37b| ) follows from the convergence at i?^ = of the integral term in 

h{B^) = 2m + B - [ d^U e-'^'-^fiU) . (4.38) 



This in turn proceeds from the behaviour of f{U) when U or f/ — > oo which may be inferred 
using ( ^4.301 ) and ( [4.27a| ) namely 



fW) .r ^r<iQ=V«^). (4.39b) 



u-^oo kU"^ u^oo kU^ Jm 

The second boundary condition follows from the fact that in the B ^ oo limit, V{B'^) vanishes, 
and f{B) tends to /o(5). Plugged into (|]33D and this results in (|]37|). 
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The spectrum. Going to -B-space in the expression for the radiation density ( 4.28|) , we obtain 

K=0 



dl 



2a 



UJ 



duj dz TT^A 



Re 



(2^ 



UB).f{-B) 



(4.40) 



The subtraction term (k = 0) given by the i?-image of the trivial solution ( [4.31 ) corresponds to 
h{B'^) = 0. Bearing this in mind, we substitute ( [4.33| ) and make use of ( |4.36| ) to arrive at 



dl 



CO 



dw dz 



2a 

— Re 

Att 



dB^ 1 - V{B^ 



2m 



52 



KB') 



Now, invoking the differential equation ([4.37]) leads to a surprisingly simple result: 



dl 



Aa 



Re{/t /i'(0)} . 



(4.41) 



(4.42) 



duj dz Att^ 

Thus the determination of the spectrum is equivalent to the following mathematical problem: 
calculate /i'(0) for h the solution of ([4.37[). 



4.3 Solution in the small k limit 



The previous discussion is valid for general V^s and therefore the differential equation ( [4.37[ ) is 
easily treated in the limit k <^ 1 by applying the WKB method. The (appropriately normalized) 
WKB-solution reads 



h{B' 



2m 



-V'{B^)B' 
L - V(52) 



exp < —] 



Ak Jo 



dB 



12, 



1 - V{B 



I2\ 



] 



(4.43) 



As long as (1 — V{B'^))/B'^ ^ —V'lB"^) has a finite B^ —>■ limit, this solution can be applied 
down to B^ = 0. This corresponds to the case when the mean squared transverse momentum 
{qj_) is well defined: 



d^Q V{Q^) = 1 ; 



-4 



d 



dB^ 



d'Q V{Q' 



[iBQf 
2! 



d^QQ^ViQ-") = 



Ml 



< 00 



Evaluating /i'(0) and substituting into ( [4.42| ) yields the radiation density 

dl ' 



uo- 



^^2K{-AV'm = -1 

ATT 71 



(4.44a) 
(4.44b) 

du dz Att ^ ^ " vr V XE'^ ' " ' ^ ^ 

which coincides with the original Migdal's result P,^. In Appendix ^ we make a closer contact 
with the derivation of Migdal, based on the random walk picture. 

In the case of interest here of Coulomb interactions, the WKB approximation ( 4.43| ) is not 
suitable, since V^B"^) developes the logarithmic singularity near B^ = 0. This case is dealt with 
in Appendix B, the result of which is corroborated by a more rigorous calculation of h'{0) due 
to Chadan, Martin and Stubbe |T^. The final answer reads 



K < 1 



dl 



a 



LU 



duo dz 



— In- 

Att V K 



K < 1 



(4.46) 
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5 Conclusion 



The aim of this paper has been to revisit the LPM effect in QED. The suppression of the radiation 
spectrum due to destructive interferences between radiation amphtudes induced by multiple 
scattering on static Coulomb centres has been studied. In the soft photon approximation and in 
the limit of large electron energy and infinite medium, the radiation spectrum has been shown 

to depend on the single parameter k = -^-^ which characterizes the coherent nature of the 
effect: the LPM suppression appears when k < 1. 

The present approach corrects the derivation given in |Q for QED. The result for the spectrum 
is only slightly different from within a minor change in the logarithmic factor. 

The dominant interference terms are due to centres separated by a distance of the order of the 

coherence length uX = —^X much larger than the mean free path A, which corresponds to large 

formation times of the radiated photon. The spectrum is determined by photon angles which 

become as large as = U'^-^t: ~ —;=-=^, compared to a typical scattering angle of order — . The 

Jk E^ 



fact that f/^ becomes as large as means that the region B'^ ~ ^/k <C 1 is the dominant region 
in "impact parameter" space. In case the Fourier transform of the normalized scattering cross 
section V^(i?^) has a finite derivative V^'(O), corresponding to a scattering potential decreasing 



faster than at large Q, the radiation spectrum is proportional to y (o'l) = y ~4y(0)yu2. For 

a potential which is Coulombic at short distances V'{B^\ ~ — In — and the logarithm gives 

4 

an enhancement factor Jin —= in the radiation spectrum. For a general potential ( ^4.42| ) gives 



the radiation spectrum in terms of the solution of a "Schrodinger" equation ( [4.3?| ). 

Finally, the present procedure can be naturally generalized to high temperature QCD in the 
spirit of [ll]. The result of [Q] for QCD should be modified accordingly, which gives for the 
radiated gluon spectrum 

^ , , = -p. T-W'^QCD In , (5.1) 

duj dz 271 Xg V kqcd 

where Cr = Cp {N^) for a fast propagating quark (gluon). This result is valid for kqcd = 
Xgjj? /uj <^ 1, where is the mean free path of the gluon. 



The derivation given in the present work may be extended to finite length media [|13[, which 



is important for future phenomenology. Last, since the method of derivation of the radiative 
energy loss induced by multiple scattering is independent of the detailed form of the scattering 
potential, provided it satisfies some requirements at large momentum transfer, a generalization 
to cold nuclear matter seems to be possible 0,113 . 
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A Gaussian interaction in the limit <^ 1 



For comparison it is instructive to relate our derivation of the radiation spectrum to the one 
originally performed by Migdal . Instead of averaging over transverse momenta and longitu- 
dinal coordinates as described in ( p.l9| ) a Gaussian probability density is used as a solution 
of an underlying Fokker-Planck equation , which reflects the assumed random walk nature 
of multiple scattering. When fixed time steps equal to A are considered the normalized density 
reads 



n 



exp(-^) 



(A.1) 



where /x^ is identified as the squared average scattering transverse momentum, /x^ = < >. 

Starting from ( [4.20|) , and averaging with ( [A.l|) , one gets ( [4 .241) where V is now the Gaussian 
interaction 

ViQ') = -expi-Q'), (A.2) 



TC 



and the function ip{U'^) is just given by the phase ipiU"^) = exp(iKf/^) ~ 1 + inU'^ for k -C 1 
which indeed coincides with ( [4.24b| ). In order to obtain the soft radiation spectrum we follow the 
steps of section noting that in i?-space V^'(O) = —1/4. The result is already stated in ( ^4.45| ). 
The absence of a logarithmic dependence on n is qualitatively understood since scattering at 
large transferred momentum is exponentially suppressed . 

The averaging prescription of a random walk ([A.l|) is also used by Blankenbecler and Drell 
IP, but expressed in terms of the transverse electric field off which the charged particle scatters. 



B The spectrum for the Coulomb case in the limit k <ti 1 

Here we derive the value of h'{0) in the Coulomb case. The Fourier transform of VlQ"^) is 



V{B' 



B K^{B) 



TT (Q2 + 1)2 

~ l-—\n^ + 0{B^) 
4 52 ^ ^ 



with the modified Bessel function Ki 15 



The differential equation (|4.37|) becomes 

4m h"{B^) + 



For <^ 1, h satisfies the approximate equation 



h"{B^^ 



1 



The approximate solution 



h{B'^) ~ Cexp 



In ( J, ) KB'' 



16k 



(B.i: 



(B.2) 



(B.3) 



(B.4) 
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may be seen, by implementation in (p.3| ), to be valid in the restricted region 




— < 52 < 1 



In — > 1 . 



(B.5) 



The constant C and the sign of the exponent in ( [B.4| ) are fixed by continuity when looking at 
the asymptotic forms of the solution h of (p.2| ). For 5^ ^ 1 [jl5| 



hiB') ~ C'BH\ 



(1) 



4^2 V 



C'WB exp 



-l+i B^ 



4^2 V K 



(B.6) 



which satisfies h{oo) = 0. The constants C and C" are of order unity. 

Extending the forms ( [B.4 ) and (B?^) to their limit of validity i?^ ~ 1, it is clear that when 
K ^ 1, we have to choose the minus sign in ( p.4| ). 

For 5^ — s> 0, the solution of ( [B.3|) may be written as 



h{B^) ~ 2i7f + h'{0)B^' 



TT 



B^ 



This expansion is valid for 



16 K 



In 



1 

52 



0{B^) . 



(B.7) 



B.8) 



Assuming that K\h' 



In (fi:|/i'(0)| 

^ni^ (this will be checked a posteriori in ( |B.11| )), ( [B.7| ) is valid for 



(B.9) 



Using (|B.5|) and ( p.9|) , it appears that the forms (p.4| ) and (p.7|) must coincide in the region 




(B.IO) 



Thus we have 



C 
h'{0) 



2171 . 



I -I 



4 V 

The final result for the radiation spectrum is obtained from ([4.42|) 

dl 



uj- 



duj dz 



K<C1 



In - . 



(B.ii: 



(B.12) 



The dominant impact parameter region which determines this spectrum is given by 

^^-v^, (B.13) 

corresponding to 



1 



^ > 1 



:b.i4) 
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